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Abstract 

The precise investigation of the W~^W~ production at the e~^e~ International 
Linear Collider(ILC) is of crucial importance in probing the couplings between massive 
vector gauge bosons and discovering the signature of new physics beyond the standard 
model(SM). We study the full one- loop EW effects on the observables, such as, the 
total cross section, the differential cross section of the invariant mass of Vl^-pair, the 
distribution of the angle between Vl^-pair, the production angle distributions of W- and 
.Z^*^-boson, the distributions of the transverse momenta of final W- and Z'^-boson, and 
the forward-backward charge asymmetry of Ty~-boson. Our numerical results show 
that the EW relative correction to the total cross sections((5e'uj) varies from —17.6% to 
—5.3% when mn = 120 GeV and ^/s goes up from 300 GeV to 1 TeV. 



Keywords: W'^W production, W-boson pair, electroweak radiative 
corrections 

PACS: 12.15.Lk, 12.38.Bx, 14.70.Hp, 14.70.Pw 



1 



I. Introduction 



The Higgs mechanism plays an important role in the Standard Model (SM) [I112]- It describes 
that the longitudinally polarized components of the physical Z^- and W^^-bosons eat the 
hidden degrees of freedom of the Higgs field. The SU{2) x f/(l) gauge invariance provides 
stringent constraints on the strengthes of triple and quartic gauge couphngs. The accurate 
measurements of these couplings could provide the information about the electroweak (EW) 
symmetry breaking. 

The multiple gauge boson productions are suitable for probing the self-coupling properties 
of the gauge bosons, and would give a crucial test of the non-Abelian structure of the SM. 
If the measured cross section is in agreement with the SM prediction, we can put a severe 
constraint on new physics. On the contrary, if there really exist gauge boson anomalous 
couplings, it would generally lead to sizable effects on the EW observables. Therefore, 
probing gauge couplings and searching for possible anomalous contributions due to the effects 
of new physics is one of the most important tasks of the present and future high energy 
experiments. 

Among all the gauge boson self-couplings, the triple gauge couplings(TGCs) of the neutral 
EW bosons Z^, 7 and the charged bosons have been well measured at the LEP2 The 
e^e^ — W~^W~ process at the LEP2 was measured not only for determining the W mass, 
but also for probing the charged TGCs[4J. To match the experimental accuracy, the one-loop 
level EW corrections to e+e" W~^W~ and e+e~ W'^*W^* 4/ were calculated in 
Refs. 13, E] . The logarithmically enhanced two-loop electroweak radiative corrections to the 
differential cross section for ly-pair production at the ILC up to the second power of the large 
logarithm were also provided in Ref.[7]. The experiments at LEP2 demonstrated that the 
SM expectations are in good agreement with the experimental data within a few percent 
If the colliding energy is larger than the threshold of Z'^-boson pair production, the Z°-pair 
production process can be used to probe the neutral TGCs. At the Fermilab Tevatron, the 
CDF and DO collaborations performed also some experiments about the diboson production 
in pp collisions at ^/s = 1.96 TeV, and presented the limitations on anomalous TGCs in 
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Ref.[H]. 

Triple massive gauge boson production processes, such as Z^Z^Z^ and W~^W^Z^ produc- 
tions, will be investigated at the Large Hadron Collider(LHC) and future International Linear 
Colhder(ILC). These processes can be used to probe the quartic gauge couplings(QGCs). In 
Ref.[n], the precise predictions for the VVV productions at hadron colliders were provided. 
It shows that the QCD corrections increase the Z^Z^Z^ and the W^W~ Z^ cross sections 
at the LHC by about 50% and 70%, respectively. Therefore, any quantitative measurement 
of the concerned gauge couplings at hadron colliders will have to take QCD corrections into 
account. 

Compared to hadron machine, e^e~ linear collider has the advantage in performing exper- 
imental measurement with a particularly clean environment. Actually, our present knowledge 
about particle physics came from both types of colliders. For example, the Z° and mas- 
sive gauge bosons were firstly discovered at a hadron collider, but their detailed properties 
and roles in the SM theory were from the LEP experiments. Therefore, lepton and hadron 
colliders are always complementary machines. 

The future ILC is an efficient machine for precise experiments with e^e~ colliding energy 
range of 200 GeV < y/s < 500 GeV in the near future. It would be upgraded to y/s ~ 
1 Te\^[TD]. This machine has sufficient energy to produce multiple massive vector bosons, 
and would be ideally suited to precision studies of the self-couplings of the vector gauge 
bosons. For example, the reactions e'^e" — ?■ Z^Z^Z^ and e+e^ — t- W^W~Z^ are very 
important processes at the ILC for probing the quartic massive gauge couplings with high 
precision. The process e"^e^ — > Z^Z^Z^ can be used to provide some informations about 
the anomalous Z^Z^Z^Z^ coupling, and its one-loop EW corrections have been calculated 
in Ref.[TTj. The phenomenology of the process e+e~ — )■ W^W~ Z^ at the leading order (LO) 
was studied in Ref. [12] . In order to match the experimental accuracy, it is necessary to take 
into account the EW radiative corrections in the theoretical predictions. 

In this work we calculate the complete one-loop EW corrections to the process e"'"e" — t- 
W^W~ Z^ in the SM. The paper is organized as follows: In Section 2 we describe the calcu- 
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lations of the leading-order (LO) cross section and the full 0{aew) EW radiative corrections 
to the e^e^ — t- W^W" process. In Section 3 we present some numerical results and 
discussion. Section 4 summarizes the conclusions. 

II. Calculations 

We adopt the 't Hooft-Feynman gauge in the LO and next-to-leading order(NLO) calcula- 
tions, except when we verify the gauge invariance at the LO. The FeynArts3.3 package[13j 
is employed to generate the Feynman diagrams and their corresponding amplitudes. The 
reductions of the amplitude are mainly implemented by using FormCalc5.3 programs |14]. 
Since the contribution from the Feynman diagrams involving — e"*" — or — e+ — e~ 
coupling is negligible due to the Yukawa coupling strength being proportional to the related 
fermion mass, we do not involve these graphs in our calculation. Then there are twenty 
Feynman diagrams for the process e^e~ — )■ W^W^ at the tree-level (shown in FigJT]). We 
denote the process e^e~ — )■ W^W~ Z^ as 

e+(pi) + e-(p2) H^+(P3) + W'{p,) + Z°(p5), (2.1) 

The differential cross section for the process e^e~ W^W~ Z^ at the LO is then obtained 
as 

ddtree = \Mtree?d^^, (2.2) 

spin 

where M.tree is the amplitude of all the tree-level diagrams, and the factor ^ is from taking 
average over the spins of the initial particles. The three-particle phase space element ^$3 is 
defined as 



In the EW NLO calculation we take the definitions of one-loop integral functions as 
presented in Ref. [15]. The complete EW one-loop Feynman diagrams include 3510 graphs, 
and we organize them into self-energy (1280), triangle (135 7), box(605), pentagon(140) and 
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counterterm(128) diagram groups. Some of the pentagon graphs are depicted in Figj2]as a 
representative selection. We adopt the dimensional regularization(DR) scheme [12] to regu- 
larize all the soft IR and UV divergencies, where the dimensions of spinor and space-time 
manifolds are extended to D = A — 2e, to isolate the UV and IR divergences. The collinear 
IR singularities are regularized by keeping finite electron/positron mass. The Cabibbo- 
Kobayashi-Maskawa(CKM) matrix is assumed to be identity matrix in our calculation. We 
adopt the definitions for the relevant renormalization constants as presented in Ref . [12] . Us- 
ing the on-mass-shell conditions [T7], the relevant renormalized constants can be expressed 
as [15] 

5Zf = -i?eS^(«)K) - ml^^Re [Sf (p^) + Sf (p^) + 2Sf (/)] |,._., 
dmy^, = Rellj, [m^), bZw± = -Re — — — \p^=m^^, oZaa = tt^ — |p2=o, 



tie —7 — La^mi, d^zA — ^ — , c)Z^2 — — 2/te ^ 



For the derived charge renormalization constant and the counterterm of the parameter s^y, 
we have[ 



5Z, = --6Zaa-~5Zza. — = -^T^e p'^'^7^^ - ^'^' Y"' ] ■ (2-5) 
2 cw2 sw 2s^ V "^VF J 

The reductions of the vector and tensor integrals are done exactly by using the approach 

presented in Refs.[T8l [T9|. The numerical calculations of the scalar one-, two-, three-, four- 

and five-point integral functions are processed according to the expressions presented in 

Refs.[Tni 1201 EI]- The calculations are carried out by using LoopTools-2.4 package [H] [22] 

and our independently developed programs for the calculations of scalar, vector and tensor 

five-point integrals with the approach presented in Ref. [IH] separately, in order to cross check 

for possible numerical instabilities. The virtual contribution of 0{a\^) to e'^e~ — )■ W^W'^ Z^ 

process can be expressed as[ 



A(T,ir = atree^vir = ^pj/l / ^^$3 J] Re (^A^tree-M^) , (2.6) 

' spin 



where pi is the c.m.s. spatial momentum of the incoming positron. A^vir represents the 



renormahzed amphtude of one-loop Feynman diagrams. 




According to the Kinoshita-Lee-Nauenberg (KLN) theorem |24]. we should consider the 
contribution of the real photon emission process e^e^ — i- W^W~ Z^'-y in order to get the IR 
safe observables for the process e^e~ — )■ W^W~ at the NLO. There includes 148 tree-level 
Fynman diagrams for the photon emission process e^e~ — ?• W~^W~ Z^'j . In the calculation 
of this process, we adopt the phase-space-slicing (PSS) method [25] to isolate the soft photon 
emission singularity. We divide the photon phase space into two parts: If < 6sEb, it's 
called soft photon region. If > 6sEh, it's in hard photon region. Then the cross section 
of the e~^e~ — )• W'^W" Z^'y process can be expressed as 

^O'real = ^(^soft + ^(^hard = O'treei^soft + ^hard), (2.7) 

where only the term Aagoft includes soft IR singularity. Theoretically, both Aasoft and 
AcThard should depend on the arbitrary soft cutoff 6s, but the total EW one-loop correction(Acrto() 
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and AcTr-eaZ should be cutoff 5s independent. 

In dealing with the soft IR divergencies, we introduce a fictitious small photon mass(m^) 
for the internal photon lines of loop diagrams, and reproduce the soft IR divergent integrals 
upon the replacements of 

7r-^-r(e,^) ^ ln(m2), ^.-^m'^hll ^ i ln\m'^), (2.8) 

where is chosen with a sufficiently small value, but not too small to induce numerical 
instabilities. After doing the replacements of (12. 8p for the IR divergent integrals, we give up 
the use of DR scheme for the case of IR divergences and adopt the massive photon scheme 
with a fictitious photon mass as regulator. That replacements are also done in treating with 
the soft IR singularity for the process e~^e~ — )■ W~^W~ Z^j before integrating over the phase 
space for the emitted photon. Generally we take a small value for 6s in the calculation for 
the process e~^e~ — )■ W~^W~ Z^'j . The terms of order 6s in Aasoft can be neglected and 
the Aasoft can be evaluated analytically by fixing a small photon mass value. In the hard 
photon phase space region, Aahard is calculated with photon mass being set to zero. After 
regularizing the soft IR divergencies with massive photon scheme, the UV finiteness of the 
whole contributions from the virtual one-loop diagrams and counterterms has been checked 
both analytically and numerically in DR scheme. 

If the IR singularity in the soft photon emission process is really cancelled by the virtual 
photonic corrections, the independence of Acrtot(= AcTmrtuai + AcTreai) on the cutoff 6s and 
fictitious photon mass m^, should be demonstrated in numerical calculation. The phase 
space integration for hard photon emission process e^e~ — W~^W~ Z^'-f can be computed 
directly by using the Monte Carlo method, because it is UV and IR finite. In practice we 
perform this integration in hard photonic region by using our in-house 2 — 4 integration 
program based on Monte Carlo integrator Vegas. Finally, the EW NLO corrected total cross 
sect ion ((Ttot) up to the order of 0{a^^) for the e+e^ — )■ W^W~Z^ process is obtained by 
summing the 0{al^) Born cross section(cr(ree), the 0{aj^) virtual cross section(AcrOTr), and 
the 0{a^^) cross section of the real photon emission process e"^e~ — )■ W~^W~ Z'^-y (Aareai)- 

(1 + 5e«,) , (2.9) 

8 



entree (/&)(FeynArts) 


octree (/&)(FeynArts) 


Crtreeifb) (CoiIipHEP) 


atreeifb) (CompHEP) 


Feynman gauge 


Unitary gauge 


Feynman gauge 


Unitary gauge 


35.810(4) 


35.810(4) 


35.80(2) 


35.80(2) 



Table 1: The comparison of the numerical results of the LO cross sections atree 
for e~^e~ — )■ W^W~ process with y/s = 500 GeV in Feynman gauge and uni- 
tary gauge, by taking the related input parameters as in Eqs. fl3.ip and using Fey- 
nArts3.3/FormCalc5.3 and CompHEP-4.4p3 packages separately. 

where Sew is the full 0{aew) EW relative correction. 

III. Numerical results and discussion 

In the following we perform the numerical evaluations at the LO and EW NLO in the a^w- 
scheme and the relevant input parameters are taken as |23j : 

mz = 91.1876 GeV, mw = 80.398 GeV, sin^ Qw = I - = 0.222646, 
mu = md = 66 MeV, = 104 MeV, = 1.27 GeV, ^ /g ^\ 

mfe = 4.2 GeV, mt = 171.2 GeV, rUe = 0.510998910 keV, 
= 105.658389 MeV, = 1776.84 MeV, 

where we use the effective values of the light quark masses [rriu and rrid) which can reproduce 
the hadron contribution to the shift in the fine structure constant ae«,(m|)[26j. In the LO 
and NLO calculations we take the fine structure constant ae«,(0) = 1/137.036 as input 
parameter. 

For the numerical verification for the correctness of our LO calculation, we use both 
CompHEP-4.4p3 and FeynArts3.3/FormCalc5.3 packages to calculate the LO cross section 
of process e"^e~ — W^W~ by adopting 't Hooft-Feynman gauge and unitary gauge sep- 
arately. The numerical results are listed in Table [H It shows they are in good agreement. 

As we mentioned in above section if our NLO calculation is correct and the IR divergency 
is really cancelled, the total cross section should be independent of and 5s- In fact, our 
calculation shows when the fictitious photon mass varies from 10^^^ GeV to 10~^ GeV 
in conditions of 5s = 10^^, mn = 120 GeV and ^/s = 500 GeV, the numerical results for the 
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Figure 3: (a) The 0{a'^^) contribution parts of cross section for e~^e~ — )■ W^W~ process 
as the functions of the soft cutoff Sg in conditions of = 10~^ GeV , uih = 120 GeV and 
^/s = 500 GeV. (b) The amphfied curve for Aatot of Fig.3(a) versus 6s including calculation 
errors. 

cross section correction Aatot = Aareai + Aa^r, are in mutual agreement up to ten effective 
digits. The independence of the total EW NLO contribution to e~^e~ — )■ W~^W^Z^ process on 
soft cutoff 6s is demonstrated in Figs|31^a,b), where we take = 10~^ GeV, rriH = 120 GeV 
and = 500 GeV. The amplified curve for Aatot in Fig|31^a) is depicted in Fig|21^b) 
including calculation errors. FigsJ3]^a,b) show that although both Aa^ij- + Adsoft and Adhard 
are strongly related to soft cutoff 6s, the total EW NLO contribution Aatot = Aa^ij. + Aerial 
is independent of the cutoff 6s within the range of calculation errors as expected. In further 
calculations, we fix = 10^^ GeV and 6s = 10~^. 

In FigJD^a) we depict the curves for the LO and EW NLO corrected cross sections as the 
functions of colliding energy ^/s with mfj = 120 GeV. FigJD^b) shows the corresponding 
relative corrections (5e«) = ^oi the data drawn in FigJD^a). We find from FigsJH^a,b) 

that the LO and EW NLO corrected cross sections are sensitive to the colliding energy y/s 
in the range of y/s < 800 GeV, and the LO cross sections are suppressed by the EW NLO 
corrections in the whole ^/s range plotted in FigJU^a). FigJD^b) shows that the absolute 
relative correction can be very large in the vicinity where ^/s approaches to the threshold 
of W+W^Z° production. That effect comes from the Coulomb singularity in the Feynman 
graphs involving the instantaneous virtual photon exchange in loop which has a small spatial 
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Figure 4: (a) The LO and EW NLO corrected cross sections for the process e^e~ — ?■ 
W^W~ as the functions of colhding energy ^/s with uih = 120 GeV . (b) The corre- 
sponding EW relative radiative corrections versus a/s. 

momentum. To show the numerical results more exactly, we list some representative numer- 
ical results of the LO, EW NLO corrected cross sections(o"(,.ee, c^ot), the EW NLO correction 
to the cross section {/S.atot = c^tot — o'tree) and the EW relative correct ion (^e^^ = Aatot/crtree) 
in Tabled There we take tuh = 120 GeV, 150 GeV and ^ = 300 GeV, 500 GeV, 800 GeV, 
1000 GeV, separately. The results in Table |2] show that both the LO and NLO corrected 
cross sections for the e~^e~ — )■ W~^W~Z^ process are insensitive to Higgs boson mass. From 
FigJD^b) and Table [2] we can see when mn = 120 GeV and -y/i goes up from 300 GeV to 
1 TeV, the EW relative radiative correction 6ew varies from —17.6% to —5.3%. 

Because the distribution of transverse momenta of is the same as that of in 
the CP-conserving SM, we provide only the distributions of at the LO and EW NLO 
in Figjni^a). The differential cross sections of transverse momentum of Z°-boson at the LO 
and up to NLO (^r^ and ) are drawn in Figin](b). In these two figures we take mn = 
120 GeV and ^/s = 500 GeV. We can see from FigsJS^a,b) that the EW NLO corrections 
generally suppress the LO differential cross sections especially when (pf ) > 100 GeV. 

We take the orientation of incoming electron as the z-axis. The 6w- (or 6z) is defined 
as the "-boson (or Z°-boson) production angle with respect to the z-axis. In Figs|ni^a,b) 
we present the LO and EW NLO distributions of cosines of the pole angles of W~- and 
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(GeV) 


mniGeV) 


<^tree{fb) 








300 


120 


2.9457(2) 


2.427(2) 


-0.519(2) 


-17.62(7) 


300 


150 


3.1605(2) 


2.633(2) 


-0.527(2) 


-16.67(6) 


500 


120 


35.810(4) 


33.51(5) 


-2.30(5) 


-6.4(1) 


500 


150 


36.035(4) 


33.85(5) 


-2.19(5) 


-6.1(1) 


800 


120 


52.34(1) 


49.70(6) 


-2.64(5) 


-5.0(1) 


800 


150 


52.46(1) 


50.10(6) 


-2.36(6) 


-4.5(1) 


1000 


120 


53.21(1) 


50.37(7) 


-2.84(7) 


-5.3(1) 


1000 


150 


53.28(1) 


50.70(7) 


-2.58(7) 


-4.8(1) 



Table 2: The numerical results of the LO, EW NLO corrected cross sections, the 
EW NLO correction to the cross section (Aatot) and the EW relative correction 
(Sew) for the process e+e" — )■ W^W~Z^ , in conditions of mn = 120 GeV, 150 GeV 
and x/s = 300 GeV, 500 GeV, 800 GeV, 1000 GeV, separately. 




Figure 5: Distributions of the transverse momenta oi W - and Z°-boson for the e'^e' 
W+W'Z^ process at the LO and EW NLO with = 500 GeV and tuh = 120 GeV. 
for l^~-boson, (b) for Z°-boson. 
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Figure 6: Distributions of the cosine of the W -boson (Z^'-boson) production angle with 
respect to z-axis for the e~^e~ — )■ W^W~Z^ process at the LO and EW NLO with ^/s = 
500 GeV and tuh = 120 GeV. (a) for '^p^^^, (b) for 

boson(cos 9]y- and cos 9z) respectively, in conditions of y/s = 500 GeV and ttih = 120 GeV. 
Both LO and NLO curves in Figj6](a) show that the produced VT'-boson declines to go out 
in the forward hemisphere, while Fig|6]^b) demonstrates that the LO and NLO distributions 
of the outgoing Z°-boson are symmetry in the forward and backward hemisphere regions. 

The distributions of the M^-pair invariant mass M\y]Y at the LO and EW NLO are shown 
in FigI7](a), and the differential cross sections of the cosine of the angle between the produced 
1^-pair at the LO and EW NLO are presented in Fig|7](b) where we take rriH = 120 GeV 
and ^/s = 500 GeV. We can see from the FigJTl^a) that there is an enhancement in the 
relatively large M^-^^ region(from 350 GeV to 400 GeV) for each of the LO and NLO 
distributions, and the EW NLO correction suppresses significantly the LO differential cross 
section dcrio/dMyyiY in this region. FigjTl^b) shows the LO and EW NLO distributions of 
cosine of the angle between the produced VT-pair. And we can see from the figure that the 
produced 1^-pair prefer to go out almost back to back, that leads to the Myy^ having the 
tendency to distribute in large value region. That is why we see an enhancement in the large 
MwiY region for each of the LO and NLO distributions as shown in Fig|71^a). 

On the analogy of the definitions of the LO and NLO forward-backward charge asymme- 
try for top-quark in Ref. |2^ , we define the LO and NLO forward-backward charge asymme- 
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Figure 7: (a) Distributions of the invariant mass of VT-pair at the LO and EW NLO when 
rriH = 120 GeV and a/s = 500 GeV . (b) Differential cross sections of the cosine of the 
angle between the produced py-pair at the LO and EW NLO with ttih = 120 GeV and 
y/s = 500 GeV. 

tries of iy~-boson as, 

aW- _ ^tree aW' _ ^tree . ^^tot ^^tot 

^FB,LO — + ' ^FB,NLO ~ + \ - + 

^tree ^tree \ ^tree ^tree 

The explicit expressions for a"^^^ are defined as 

(^tree = (^treeiVw- > 0) ± atree{yw~ < 0), (3.3) 

where i/w- is the rapidity of "-boson, the notations (TtreeiVw- > 0) and crtreeiuw- < 
0) represent the cross sections for the produced iy~-bosons in the forward and backward 
hemispheres at the LO respectively. The forward direction is along the orientation of z- 
axis. ^(J^t denote the EW NLO corrections to the cross sections cr^.f.e- Iii the conditions of 
^ = 500 GeV and uih = 120 GeV, we get A^^^^ = 50.98(2)% and A^j^^^ = 53.28(3)%. 
Both the LO and NLO results show that most of the "-bosons are produced in the forward 
hemisphere, that feature has been already demonstrated in Figl^a). 

IV. Summary 

The W^W~ production via electron-positron collision at the ILC is an important process 
not only in probing the non-Abelian structures of the SM, but also in finding new physics. 
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In this report we have shown that the phenomenological effects due to the one-loop EW 
radiative corrections, can be demonstrated in the e^e~ — W^W' process for all colliding 
energies ranging from 300 GeV to 1 TeV at the ILC. Our results show the EW one-loop 
radiative corrections significantly suppress the LO cross sections, and the relative correction 
to the cross section varies from —17.6% to —5.3% when itlh = 120 GeV and ^/s goes up 
from 300 GeV to 1 TeV . We can see the obvious effects of the EW NLO correction on the 
physical observables, such as, the distributions of the transverse momenta of final Z^- and 
ly-bosons, the differential cross section of the invariant mass of VT-pair, the distribution of 
the angle between ly-pair, the production pole angle distributions of W~- and Z'^-boson, 
and the forward-backward charge asymmetry of iy~-boson. 
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